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In an earlier paper we reported that an empirical formula of the electrical conductivity a — 
e^natomTo/(?nG) agrees with experiments within ~ 20% for the most of pure elemental metals at 
room temperature ranges. Here tq — h/ksT is assumed for 'all' metals and G is summed electron 
numbers in each atomic shell: e.g. G—6 for Cr(3d^4s^). In this paper, we find that the above to 
can be deduced if the autocorrelation time of electron fluctuating velocity in a simple fluctuation 
dissipation theorem is converted to 2 <AE><At> /ksT, and if this < AEx At> is assumed 
equal to h/2 of the Heisenberg's minimum uncertainty. This corresponds to the cloesest approach, 
or head-on collisions. Independent from this, we find that Tidcctron = natom is appropriate for a in 
most elemental metals . In discussing temperature dependence of a ~ T~^, besides use of the Debye 
temperature unit (©d), the temperature unit of 0p = fiUpi/kB is found to be equally acceptable. 
Here LUpi is the ion plasma frequency, only depending on ndcctron unlike somewhat ambiguous 0o- 

PACS numbers: 72.15Eb, 72.15Lh, 72.10.Bg, 65.40.Ba 
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1. Introductionn 

In the previous paper [3 (Paper I), we found empiri- 
cally that absolute experimental values of the electrical 
conductivity a for most of the pure elemental metals at 
room temperature can be well reproduced (~ 20%) by 



^ = T^'^O, 

niG 

where the relaxation time tq is assumed to be 

h 



To 



(1) 



(2) 



for 'all' metals. Here riatom is the number density of atom 
in each metal and m is the true electron mass. The factor 
G is integer and is summed numbers of outer electrons in 
the electron configuration of each atom such as G = 1 for 
Cu(4si), G=l+2=3 for ln'^^{5s^4p^), or G=l+4=5 for 
Nb"'^(3d^4s^). These G values are empirically selected 
to match the observed a. Note that we treat metals of 
normal geometrical size at ordinary pressure. 

In contrast to our formulae, the conventional Drude 
form a = ne^r/m* requires the values of Z = n/riatom 
for n (the electron number density), m* (the effective 
electron mass) , and in particular r (the relaxation time) . 
This is always presented in well-read standard textbooks 
0-11], but unfortunately they do not present conclu- 
sive ways how to calculate 'absolute values' of cr. As 
noted in paper I, though Lifshitz and PitaevskiiQ and 
AbrikosovQ suggest the form of tq ~ h/k^T, how ac- 
curate their '~' signs is not clear, nor they treat non- 
monovalent metals. These issues are of fundamental im- 
portance, but they are forgotten or ignored as research 
works in the world of solid state physics today. 

Now, examples of the ratio of tq from eq.(2) to Tobs 
determined from eq.(l) using the experimental resistiv- 
ity and G = 1 are ro/Tobs=0.98(Na), 0.97(K), 0.99(Rb), 



1.05(Cu), and 0.95(Au) at room temperature ranges 
(hereafter Troom)- Because of this remarkable closeness 
to unity, it is tempting to assume that there may be an 
extremely simple explanation for tq — h/k^T, which is 
the major theme of this paper. 

For this purpose, first we search for the connection of 
To to the fluctuation dissipation theorem, FDT0, In 
order to reproduce To = h/ksT, we find that by con- 
verting iJ = i < mv{Q)v{t) > dt used in a simple 
FDT to H =<(A£;)2>i/2<(At)2>i/2^ this H must take 
the minimum value of h/2 expected from the uncertainty 
principle. 

Secondly, we find that the electron number density n 
which does not appear in eq.([T]) is satisfied by ti = Matom, 
i.e. Z = \, which is of course only applicable to 'conduc- 
tivity of elemental metals'. 

Since eqs.(l) and (2) reproduce the experiments quite 
well for the most metals at Troonu we can reproduce the 
observed 'absolute values' of p = l/cr at low temperature 
ranges where p cx is observed, if we employ the 
Griineisen-Bloch formula. The Debye temperature 0d 
itself used here as a temperature unit is however as is 
well-known dependent upon the temperature where it 
is determined or upon methods of derivation, and thus 
not unique. Then we attempt to use the ion-plasma 
temperature as an alternative unit, determined from 
/cb6'p = hbjpi. Here Wpj = (ne^Z^/eoA/)^/^ is the ion 
plasma angular frequency where M is the ion mass. We 
find that use of the ^p-unit for the electrical resistivity 
and specific heat Gv is equally acceptable, or even better 
in that it requires only the knowledge of n, hence unique. 



2. Interpreting RelsLxation Time 

The classical representation of the fluctuation dissipation 
theorem (FDT) for the isotropic electrical conductivity 
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a is as eq.(9.1) in the seminal paper of Kubo @ 

a — tic^tk/to, (3) 



where 



TK 



m 
1^ 



< v{0)v{t) > dt. 



(4) 



Thus in order to obtain tk 
assumed in eq.®), we need 



To (G = 1 is impHcitly 



H 



1 



< mviSl)v{t) > dt^ - (5) 



Here v{t) is the fluctuating electron velocities parallel to 
the given electric field. (Thermal fluctuation is of course 
expected without electric fields, and is on the order of 
(4/(1 — /)/A^)^/^, in unit of the average energy where 
/(e) is the Fermi distribution function and N is the total 
number of electrons in a volume V. Thus the fluctuation 
is occurring around the Fermi energy E-p with the width 
« 2k'QT. See detail in Appendix A) 

Besides an obvious relation of the autocorrelation of 
< v{G)v{t) >=< v{ti)v{t + ti) > for any ti in the as- 
sumed stationary stochastic processes, we adopt H = 
<^mv{0)v{t)At> , where At is chosen to reproduce the 
value of H. This may be not unreasonable if we con- 
sider that the autocorrelation function rapidly decreases 
as exp(— t/r) for the relaxation processes. We approxi- 
mate H =< ^mv'^At >, and interpret ([5]) as expressing 

H imu^Ai >=< AEAt >= (6) 

If mv — hk and v{t)At — Ax are used, we still obtain 

H = ^ <hk{t) X v(t)At>= <kAx>= ^. (7) 

While the minimum given in the Heisenberg relation is 
<(Ai?)2>i/2<(At)2>i/2= (^AEU^iAtU, = h/2, it 

was necessary to adopt it to match tq . 

Reversing the logic, assume that H ^mv^At > 
takes the minimum uncertainty value of ?i/2 as suggested 
in eq. ([U or ([T]), and assume that eqs.Q and dH) hold, 
we obtain tq — tk and reach tq — h/k^T. Though it 
seems rather difficult to escape from this statement, we 
later in §4 discuss relation to the case of p = 1/ct oc 
and physical meaning of eq.® at the end of §4 and in 
§6 (See also Appendix B). 



3. n/natom in Multi-Band Metals 

Though eq.(l) needs only values of riatom, naturally we 
wish to know n oi Z = n/natom- This Z value is tab- 
ulated only for 20 metals in Kittel's table while the 
other quantities are tabulated for almost all metals, in- 
dicating that Z is not easy to assign. Let us recall the 
well-known derivation a = jx/E^ for a single band as 



preparation for multi-band metals: e.g. each one of five 
'd' and one 's' bands in Cr(3d^4s^). We rewrite the cur- 
rent density as 



Jx- 



,lv,fdk'^~el^if-fo)dk', 



(8) 



using mvx — hk^jdh' = 2c?k/(27r)'^and J u^^/o^k' = 0, 
while /o and / are unperturbed and perturbed Fermi dis- 
tribution function, respectively. Given the electric field 
Ex in the x direction entering the 'steady' Boltzmann 
equation 



df df 
dt ^ hdkx 



-eEx 



dfo 
hdkx 
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we obtain from eq. 



(9) 



(10) 



by eliminating (/ — /o) and using k-independent tq of 
eq.(l2]). Because of the factor dfo/dkx, not only those 
bands in which the state density does not extend be- 
yond the Fermi energy do not contribute to a such 
as 3d^° of Cu as noted in Paper I (see e.g. ref.[l| 
Fig. 7. 12, left), but also only electrons having the en- 
ergy of (1 ± 2kaT / E-p)E-p « E-p are contributing to 
a {Ep—¥erim energy=2-7eV for Z—1). Partial inte- 
gration over dkx within the primitive cell leads to the 
second equality. The final equality leads to n, where 
n = fcp/37r^ and kp = {2mEpY/'^ /h hold. This means 
that (J / [e^TQ/m) does not depend upon the state density 
of electrons eiX Ep, nor respective bands, but it is simply 
equal to n. Note that even for the multi-bands discussed 
below, the departure from dk = Ank^dk (by a multiple 
factor of 1.2-1.6 for most metals) appears to change only 
the value of Ep The rough derivation above may be 
sufficient for the discussion below. 

For multi-band metals, we adopt the total resistivity 



(11) 



where pi = l/cri as in Aschcroft and Mermin Q (see 
eqs. (13. 20-22) and discussions followed). The average 
momentum equation can be obtained from eq. ([9]) by mul- 
tiplying by mvx/n and integrated over dv, where df /dt 
leeds to the inertia term mdVx/dt below. 



mdVx/dt ~ —eEx — tt^Vx/tq ~ 0. 



(12) 



Obviously electrons are on the average decelerated (or 
'resisted') by collisions with ions or by phonon scattering 
{—mVx/To). If there are different kinds of ions such as 
impurity, we should add up T,mVx/Ti = rnVx/rtot as the 
summed resistance (as in the Matthiessen rule). We as- 
sume that the same situation occurs among the different 
bands in elemental metals which may behave indepen- 
dently or additively. This means that for the multi-bands 
metals tq should be replaced by rtot, where rtot is given 
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by l/rtot = S(l/Ti). As a result of pi = m/^ne^Ti), the 
total resistivity is p = Spj = ra/(ne^Ttot)- Further we 
assume that every is equal to tq = H/k^T^ as can be 
surmised by the derivation for the single band. Because 
G is found to be the sum of outer electron numbers in 
electron configuration, which is nothing but numbers of 
responsible bands, l/rtot = — G/tq results, lead- 

ing to a = e^T^n/ (mG). 

Since we found a = e^Toni^to^/ (mG) in eq. ([T]) empiri- 
cally from G « Gobs J it is compatible with the above only 
when Zcff = 'T./'T-atom = 1 holds. Namely the 'effective' 
valence electron number Z^s responsible for the electrical 
conductivity is 'unity' for the most of elemental metals! 
This is consistent with G=l for Cu(3d^''2s^), while e.g. 
for the cohesive energy, use of Z(Cu)=ll may be relevant 
[Tq] . It seems thus unnecessary to introduce the effective 
electron mass m* to find a and A (thermal conductivity, 
discussed in Paper I), if not for other purposes. 

Summarizing, we assume = tq for each i-band and 
p = T,pi cx I](l/Ti) = Eil/ro = G/tq. Since this 
p — raG I (ne^To) must be equal to p — TOG/(natome^To) 
from the empirical eq.(IT]), we obtain for p of all elemental 
metals 



1. 



(13) 



Thus e.g. in Al([Ne]3s^3p), Z^s is 1 for p, different 
from Zcff = 3 for the cohesive energy, for example. 
Besides, G=l for Al is compatible with the experimental 
p (Paper I). Note that our assertion of eqs.(IT]) and ^ 
holds whether Z^s = 1 or not. 



4. Two Temperature Units 

To proceed on the theoretical interpretation of tk, we ex- 
tend the resistivity to low temperature ranges, by adopt- 
ing the Griineisen-Bloch function Fgb = ^gb(T) (see 
e.g. eq(9.62) in ref.g); Fqb « T for T.ed » 1, i^GB « 
for Tred < 1, and Fgb(1) = 1- Here T^d = T/Ot^ and 
is the Debye temperature {T^ed = ^reduced)- Thus we 
extend our formulae of eqs.(l) and (2) to the following, 
assuming that the Griineisen-Bloch formula is valid. 



P ■ 



mG 1 



]_ ^ 1 Fgb{T)/T 

Ta ToFGB{Oa)/6a 



(14) 

(15) 



Here the subscript a stands either for D or P the ion- 
plasma temperature ('D' for Debye and 'P' for plasma). 
For the high temperature range where -Fgb oc T holds, 
eq. (15) gives = tq cx T^^ and we obtain p cx T. 

Figure 1 in upper right shows the temperature depen- 
dency of the observed electrical resistivity p divided by 
temperature, p/T[r2m/K], in unit of p{0d)/Ou for 19 
metals. We adpoted the same metals used in Fig. 2 
for easy rough comparison. Here the Debye temperature 
is adopted from Kittel [4] , which is empirically deter- 
mined from the lattice specific heat at constant volume 




FIG. 1: Electrical resistivity vs. temperature for 19 met- 
als. The experimental data of {p{T) /T)/{p{0a) / 0a) are 
plotted against T/O^ in black cross sign (left a^Plasma 
and right Qf=Debye) (symbolically denoted as p/T below 
respective plots). The red plus sign plots are data from 
O.STmcit < T < Tmcit, the melting point. Down-shifted plots 
denoted as Dbyi/2 use the Debye temperature fitted around 
T = Oo/2 by use of Cv adopted from ref. 12]. Three full lines 
plot the same ratios using Griineisen-Bloch function. Though 
Fe, Co and Ni (square) show p cc T^'^, there is no discernible 
difference between 0u and 0p plots. 



Gv in low temperature regions of Gy cx T^. Because 
the tabulated temperature steps for p [ll[ are too sparse, 
i.e T(K)=[78, 273, 573, 973], logarithmic interpolation is 
made to obtain the 'observed' p{0b)- Since we are inter- 
ested whether a different unit as compared to the Debye 
temperature unit can be employed, this kind of roughness 
is not problematic. A general trend of proportionality 
from (dashed line in the middle of the figure) to con- 
stant, coinciding roughly with the Fgb /Tied (full line), is 
apparent. At higher temperature of T / 0d > 2, however, 
the departure from constancy of p/T becomes apprecia- 
ble, so that we plotted in red cross signs the observed 
points in the range 0.5Tnioit < T < Tmoit, where Tmoit is 
the melting points of each metal (factor 0.5 in 0.5Tnioit 
is arbitrarily chosen). Discussions hereafter neglect these 
points. 

For reference in the lower right, we show the p/T- 
ratio using 0d from ref. IJ] determined at the tempera- 
ture near 1 /2 of the to-be-determined Debye temperature 
from the data plots of experimental Gy . 

Because of the well-known ambiguity of adopting 
the Debye temperature (empirically from specific heat, 
acoustic speeds or p itself), and its temperature depen- 
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dency due to the various fitted temperatures as shown 
above) , we in this paper attempt to use another temper- 
ature 0p expected from the ion plasma angular frequency 



Wpi as defined below: 



9v 



hupi/ke, LOpi = {ne'Z^s/eoMY'^ (16) 



Here n is the electron number density, eo is the permit- 
tivity of vacuum (l/47r in c.g.s. Gauss unit) and Af[kg] 
is mass of each atom. The only unknown parameter in 
determining Op is the ratio of n/rig^tom = ^cff, which 
we take unity for all the elemental metals from eg. ([T3| . 
Mostly Ob/ Op is 0.9~ 1.2, while in alkali, In and Tl 
it is very close to 1/2 and near 2 in Mo, Rh, W, Os 
and Ir {Ou differs among various metals more than Op). 
The dispersion relation is w = ujpik/{k^p + fc^)^/^. Here 

kTF{= [e^D{EF)/eoY/^ = {Akp/naoy/^ < n^/^) is the 
Thomas-Fermi shielding wave number, where D{Ep) is 
the electron state density at the Fermi energy, kp is its 
wave number and oq is the Bohr radius (e.g. ref.Q in 
eqs.(17.55) and (26.4)). 

Figure 1 left plots [p{T)/T]/[p{Op)/ Op] against T/Op 
as in the Debye-plotting. Both plots either by use of De- 
bye or Plasma temperature unit are found quite similar, 
indicating that to adopt Op is not out of question, ir- 
respective of theoretical reasoning. Full black lines are 
all Fqb functions, passing through the point of (1.0,1.0), 
and of course it is identical to the ones which used 
unit. Thus if we adopt eqs.(14) and (fT5|) for a=P, namely 
adopt the temperature expected from the ion-plasma fre- 
quency as in ea. (|T6l) . most of the 'absolute' values of re- 
sistivity of elemental metals can be well reproduced with- 
out introducing empirical parameters other than G in a 
temperature range somewhat below the melting points 
at ordinary pressure. Eq.(14) with a = P might be used 
for further study, because, as is well-known, not only Fe, 
Co, and Ni depart appreciably from this in Tred = 1 — 2, 
but also below a few Kelvin many metals [l3, [Hi show 
p=const-|-const X T-' including Na, Cu, and Ag, where j w 
2 (ref.[l3| lists 10 such metals). 

In order to see the theoretical temperature dependency 
of Ta or p, we show below a simplified temperature de- 
pendency of the collision time, denoted as Tg. 

1 1 /"^ rlf) 1 /9 

— = — / (l-cos6')-27rsin6l— = — sin'*(-) (17) 
Te To Jo 47r tq 2 

for < TT and Tg = tq for > tt. This is equiva- 
lent to assume W{0) — 1/tq in eq.(3.11) of ref. 6], or 
in eq.(16.32) of ref.[3]. Here 1 — cos 61 is the weight 
factor to the direction of the given electric fields and 
2TTsm0d0/4TT is the fraction of the differential solid an- 
gle. Further we simplify that — T/ Op, instead of the 
standard sin(0/2) = q/2kp where q is the phonon wave 
number and kp is that for electrons at the Fermi en- 
ergy. Blue dash-dot lines in Fig.l left and right plot 
y(ordinate) = sin'*( Tred/2)) against X = Tred/(37r/2) 
for Trod < TT, while for > 7r/2, Y=l is kept. These 
tend to p/T cx T'^ for the low temperature, coinciding 




FIG. 2: Specific heat (usual one from lattice) vs. temperature. 
Experimental values of Cp/{3R) and Cv/{3R) are plotted 
against T/Ta in black cross marks (q — P, left and a — D, 
right). Red plus signs are those at temperatures O.STmeitT < 
Tnieit as in Fig.l. Upper 5-full lines are all the same Debye 
functions. Two bottom figures (left and right) plot Cp/Cv 
for each metal used for conversion from Cp to Cy. 



with p (X , and for the high temperature we obtain 
p oc T, which are rather satisfactory. This is found to 
be the upper envelope of the observed data points near 
X(abscissa) sa 0.3. The much simplified eq. p7|) is, if any 
merit in it, to show the dependency clearly. At the 
same time X — Trod and Y — sin^(Trod x 37r/2), and 
Y = 1 for AT > 1/3 give an identical curve, which is 
in some cases convenient since the same X value as the 
Griineisen curve is used; e.g the maximum ratio between 
two curves is 1.6 at Tred=0.2. 

Figure 2 left plots the observed Cp/(3i?) from ref.[llj| 
(upper) and Cy/{3R) calculated from Cp/(3i?) (lower, 
discussed below) against T/Op, and in the right, simi- 
lar plots against T/Op, both for 19 metals. Hrcr R is 
the gas constant. The temperature step is again sparse; 
[100, 200, 298.15, 400, 600] K. If we compare the left 
Cv/3R(plsm), meaning plasma, and right Cv/3R(Dby), 
there is no large difference between the two, and the 
Cv/3R(plsm) by the use of the ion plasma temperature 
shows smaller scatter from the Debye function (full lines) . 
Thus use of the plasma temperature as the temperature 
unit, which needs no adjustments from experiments, is 
not restricted to the resistivity as in Fig.l such that 
eqs.(14) and ([15]) with a = P may be allowed to be one 
possible standard. 

Below are some remarks. The bottom curves are 
Cp/Cv from Cp(JK-imol"i)=Cv + {3ap)^TVrn/x for 
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each metal. Here qt, = ctL(T) is the hnear thermal ex- 
pansion coefficient x is the isothermal compressibil- 
ity 0, and V-in is the volume/mol(=atomic weight/metal 
density). The down-shifted plot denoted as Cy/ZR 
(Dbyi/2) as in Fig. 1 shows naturally smaller scatter than 
the plot denoted as Cv(Dby) from the Debye tempera- 
ture determined at much lower temperature. Red plus 
plots are for Q.bT^a\t<T <T',„oit as before. Full lines are 
all the same exact Debye functions, which again passes 
(l,l)-point and the limiting value is 1.05014=1/0.95225. 
Dashed line is the limiting Cy (x T^. 

If the same correction factor Cy / Cp is tentatively ap- 
plied to (p/T)ratio in Fig.l, the corrected plot becomes 
closer to yo/T=constant for T/Tq, > 1. It is clear that 
some kind of small correction to p/T is needed, because 
the electron density will decrease due to larger separation 
between atoms when temperature increases. 

The first conclusion of this section is that use of O-p 
from the ion plasma frequency as a temperature unit 
gives reasonable temperature dependence for both the re- 
sistivity and specific heat, and it is unique as compared 
to ©D, such that eqs.(14) and (15) for a=P may be more 
appropriate. Now the applicable range of tk is restricted 
to e.g. p{T)/T « p{0p)/ Op, meaning the deflection an- 
gle be > 7r/2 as seen particularly in the sin'*(0/2)-plot 
of Fig.l. The latter implies the close encounter, or a 
small impact parameter, and suggests, as the second con- 
clusion of this section, that the minimum uncertainty of 
h/2 is maintained for tk and hence tq may well be due 
to the close collisions. These arguments could have been 
presented by employing the Debye temperature unit, but 
now by the use of 6'p these become less ambiguous. 

We note a similarity between the present close colli- 
sion and the electron in the first Bohr orbital in atomic 
hydrogen, where the quantized angular momentum 
mvr = h is conserved in the circular and closest orbital. 
In the present case of the elastic collision in a hyperbolic 
orbital, the angular momentum is still conserved, and 
its minimum absolute value is mvr — hkr ~ h/2. 



5. Electron Specific Heat 

The electron specific heat, discernible only below a few 
Kelvin or so, is given by Cye = Tr'^D{EF)klT/3, where 
D{Ep) is the number of states per energy per atom. For 
the free electron model D{Ep) = 'i/{2Ep) holds, lead- 
ing to Cye = 7r^fc|r/2i?F- The ratio Cyo,obs/C'yc,theo is 
conventionally expressed by m* jm (r7i*=the thermal ef- 
fective electron mass), and m*/m amounts to 10 or more 
for transition metals if the free electron model is used for 
the theoretical Cyc,thco (see e.g. ref.fl], table 6.2). 

However if we adopt D(Ep) from the extensive calcu- 
lation of Moruzzi et a/.jlfl] and Cyo,obs from ref.P], we 
find on the average m*/m — 1.48±0.89 for 30 metals 
down to the period starting from Rb to In. If we exclude 
six metals showing m*/m > 2.0, the average becomes 
m* /m = 1.11 ± 0.29, which is quite satisfactory. 

Values of m/m* for excluded six metals are 2.6(Sc), 



2.4(V), 2.5(Mn), 4.9(Sr), 3.1(Y), and 2.4(Nb). For ex- 
ample, Sr shows a sharp drop in D{E) near i^pil^, sug- 
gesting that D{Ep) — 0.31 could have been much larger, 
and accordingly m/m* may become closer to unity. A 
similar situation is seen in Y {D{Ep) — 1.41), which has 
a steep peak of D{E) just below Ep. 

Thus we can safely say that the theoretical under- 
standing of Cye is rather satisfactory for the elemental 
metals (except perhaps for semi-metals not included in 
Moruzzi et al.), so that we need not concern it when 
discussing other physical quantities. 

6. Discussion and Conclusion 

We found that the empirically well established tq = 
h/kp,T for the room temperature range in elemental met- 
als (Paper I) can be deduced in the following way. We 
take conversion 

I <—viO)v{t)>dt=<—-At>=<AEAt> (18) 
(see Appendix A for the second equality), and assume 



then we obtain 
m 



TK 



<AEAt>^ h/2, 



< v{0)v{t) > dt 



knT 



To- 



(19) 



(20) 



Further simplified ro-derivation is to assume from the 
outset that the relaxation time tq in the steady Boltz- 
mann equation is equal to the autocorrelation time of the 
normalized fiuctuating thermal velocity tk as in eq.(Bl). 
Namely 



eE., 



dfo 
hdkr 



I_k_ 

1 



where 



v{Q)v{t) , 
TK= / <——>dt 
Iq Vt Vt 



(21) 



(22) 



This might be regarded as a basis of derivation of the 
simple classical FDT, if so wished. Then we evaluate 
Tk, assuming eq.(19) and using eq.(18). In this case tq 
entering the steady Boltzmann equation [9] should better 
be called 'correlation time' rather than 'relaxation time'. 
At the same time the once discarded 'collision time', if 
not by all people, be retained because without notion of 
'collisions', tq ~ h/k^T cannot be understood in closest 
collisions. 

For the wider temperature ranges, we can well repro- 
duce the observed 'absolute' values of electrical resistivity 
p by eqs.(14) and ([T5|). exception being T higher than a 
fraction of Tj^^it- Here the latter assumes the Griineisen- 
Bloch function, where only the relative p values can be 
given as eq.(9.62) of ref. [2|]. 
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Further, we find that riatom used in eqs.(l) and (14) is 
equal to n, i.e. Z^g = 1 for p and for the thermal con- 
ductivity too (Paper I), though only in elemental metals, 
and not necessarily applicable to other applications (a 
short summary above eq.(13)). 

If we accept Zcs—l, the temperature unit from the ion- 
plasma frequency, ksOp = huip, becomes useful, in that 
the unique temperature value only dependent upon the 
electron density n can be defined such that the observed 
p/T and can be represented quite similarly or slightly 
better than use of 0d unit (Figs. 1-2). Regarding the ex- 
citation of ion-plasmons in eq. (jl6|) of hupi — k^Op, one 
might think of equilibration with the ubiquitous black- 
body radiation inside the metal, though one needs a care- 
ful study. 

Eq.lll]), <AEAt>^ h/2, is applicable to the large an- 
gle scattering, implying close collisions with small impact 
parameters. We mentioned at the end of §4 the similarity 
of the quantized angular momentum mvr = hin the first 
Bohr orbital. Another comparison is to use a harmonic 
oscillator (one-dimensional, 1-D) where AEAt = h/2 is 
realized in the closest oscillation with a Gaussian wave 
function, and to consider collisions with the scattering 
angle oi 9 = tt, i.e. 1-D head-on collisions (zero impact 
parameters). Difference between the two is bounded or 
unbounded motion, but with the same closest approach, 
giving the same <AEAt>— h/2. 



where use is made of p = Kk, n = 'iir'^k^^dk/k-p = 
de/2£'F, and N = nV (Gn = DNde, where D is the 
conventional state density). Multiplying eq. (|Aip by 
and using <{AEY>=<{eANY> , we find 



<E> 



3/(1 - /) de 



2N Ef 



1/2 



oc 



N 



(A3) 



The relative fluctuation eq. (jA3|) is dependent upon 
sizes (cx 1/\/N) as in the Boltzmann distribution and de- 
pendent upon the width of energy spectrum de. Besides, 
it is appreciable only near e « Ep (max[/(l — /)]=l/4) 
such that below e ~ Ep — 2k-QT , f is saturated at unity so 
that there is no fluctuations in AA^ nor AE there. Thus 



\/<(A£;)2> =< mv^/2 > 



(A4) 



should hold, where mv^ /2 is the thermal energy mea- 
sured from the center value of e = Ep. Note that if we 
adopt /(I — f)de ^ de ~ kp,T, this becomes identical to 
the problem-answer at the end of §113 in ref . [l6l| . 



Appendix B: Correlation Time in FDT 



Appendix A: Fluctuation Amplitudes 

In order to flnd the amplitude of energy fluctuations, we 
use the fluctuation of electron numbers AN in a volume 
V as given in Landau and Lifshitz [l^, eq. (113.6) 



< [ANf >=Gn/x (1-/), 



(Al) 



which is derived from partial derivative of the Fermi dis- 
tribution function / = 1/ {exp[(e — //)/A;Br] + 1} with 
respect to the chemical potential p.. This p is practically 
the same as the Fermi energy Ep. Here Gn is a scalar of 
states number given by 



G. 



3 ^h' 



As emphasized by Kubo [9] in p. 580 , unless one intro- 
duces some physical models the FDT does not give useful 
answers, even though it may be correct. This is shown 
in the 'simplest example of the independently moving 
charged particles in the classical system' (p. 585), i.e. 
eq.(4). Utilizing a normalized velocity function v(t)/vT 
where vt = {kBT/mY^^smd introducing the normalized 
autocorrelation function G, eq.(4) becomes 



°° v{0)v{t) , 
Vt Vt 



C{t)dt-- 



(Bl) 



This is almost a definition of correlation time as seen in 
C{t) — G(0)exp(— ^/tcoi ). But if it is so, tq — tk might 
be used directly in the Boltzmann equation eq.®. 
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